k = [
0.001
0.0005
0.00025
];
T = [
6.326
10.683
23.324
];
E = [
1.92629
2.54268
1.76551
];
subplot(1,2,1);
plot(k,T);
hold on;
subplot(1,2,2);
plot(k,E);
hold on;
k = [
0.001
0.0005
0.00025
];
T = [
34.879
64.466
126.958
];
E = [
0.856771
0.0255064
0.0656555
];
subplot(1,2,1);
plot(k,T);
hold on;
subplot(1,2,2);
plot(k,E);
hold on;
k = [
0.001
0.0005
0.00025
];
T = [
91.176
185.788
368.47
];
E = [
0.117327
0.0710093
0.0682191
];
subplot(1,2,1);
plot(k,T);
hold on;
subplot(1,2,2);
plot(k,E);
hold on;
subplot(1,2,1);
set(gca,'XDir','reverse');
xlabel('Time Step');
ylabel('CPU Time');
subplot(1,2,2);
set(gca,'XDir','reverse');
xlabel('Time Step');
ylabel('Max-norm Error');
legend('Euler','ClassicalRK','DormandPrince');
saveas(gcf,'pictures/winner/winner.png');
